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ABSTRACT 

We report the Fermi Large Area Telescope discovery of 7-ray pulsations from the 22.7 ms pulsar A in the 
double pulsar system J0737-3039A/B. This is the first mildly recycled millisecond pulsar (MSP) detected in 
the GeV domain. The 2.7 s companion object PSR J0737-3039B is not detected in 7 rays. PSR J0737-3039A 
is a faint 7-ray emitter, so that its spectral properties are only weakly constrained; however, its measured 
efficiency is typical of other MSPs. The two peaks of the 7-ray light curve are separated by roughly half a 
rotation and are well offset from the radio and X-ray emission, suggesting that the GeV radiation originates 
in a distinct part of the magnetosphere from the other types of emission. From the modeling of the radio and 
the 7-ray emission profiles and the analysis of radio polarization data, we constrain the magnetic inclination a 
and the viewing angle ( to be close to 90°, which is consistent with independent studies of the radio emission 
from PSR J0737-3039A. A small misalignment angle between the pulsar's spin axis and the system's orbital 
axis is therefore favored, supporting the hypothesis that pulsar B was formed in a nearly symmetric supernova 
explosion as has been discussed in the literature already. 

Subject headings: gamma rays: stars - pulsars: general - pulsars: individual (PSR J0737-3039A, PSR 
J0737-3039B) 



1. INTRODUCTION 

Pulsed GeV 7-ray emission from more than 100 pul- 
sars FI has been obse rved by the Large Area Telescope (LAT; 
lAtwood et al.l 120091) aboard the Fermi satellite, launched in 
2008 June (see The Second Fermi Large Area Cata log of 
Gamma-Ray Pulsars; iFermi LAT Collaboration! I2013L here- 
after 2PC). The current population of 7-ray pulsars includes 
objects known from independent radio or X-ray observa- 
tions and detected in 7 rays by foldin g the Fermi LAT 
data at the known ro t ational periods (e.g. .lAbdo et al.ll2009al ; 
iNoutsos et alj 120111: lEspinoza et alj 120131) . and ^-rav pul- 
sars found through direct blind searc hes of the LAT data 
(lAbdo et alJl2009bt iPletsch et al]l2012l) or r adio searches for 
pulsars in unassociated 7-ray sources (see iRay et all 120 121 
and references therein). These pulsars are energetic (spin- 
down luminosities E = 4ir 2 IP/P 3 > 10 33 erg s" 1 , where P is 
the spin period, P is its first time derivative, and / denotes the 
moment of inertia, assumed to be 10 45 g cm 2 in this work) 
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and are typically nearby. They therefore have large values of 
the "spin-down flux" E/d 2 and of the heuristic detectability 

metric V~E/d 2 (lAbdo et al.ll2010ah . 

The doub l e pu ls ar system PSR J0737-3039A/B 
dBurgav et alj 120031: iLvne et all 120041) consists of two 
radio-emitting neutron stars in a tight 2.4 hr orbit. Radio 
timing observations of the tw o pulsars provide hig h-precision 
tests of strong field gravity (iKramer et al.l 120061) . The 2.7 
s pulsar J0737-3039B has a low spin-down luminosity 
E ~ 1.7 x 10 30 erg s , three orders of magnitude smaller 
than that of the least energetic 7-ray pulsar currently known, 
making this pulsar unlikely to be detectable by the LAT. On 
the other hand, the higher £ of ~ 5.94 x 10 33 erg s" 1 — for 
which the Shklovskii effect (IShklovskiil 1 19701) is negligible 
due to the system's low transverse velocity — and the modest 
parallax distance of 1150+jgQ pc measu red with very long 
baseline interferometry observations dDeller et all 120091) 
make the 22.7 ms pulsar J0737-3039A a credible candidate 
for a detection in 7 rays. 

Searching for high-energy pulsations from PSRs 
J0737-3039A and B allows exploration of regions of 
the P— P diagram that are devoid of 7-ray pulsars (lAbdo et alJ 
l2010al) : mildly recycled millisecond pulsars (MSPs) for the 
former (mildly recycled pulsars being pulsars that were only 
partially spun-up by accretion of matter from a companion 
star), and slowly-rotating normal pulsars for the latter. 
Increasing the variety of 7-ray pulsar types helps understand 
the phenomenology of high-energy emission from these stars, 
and how emission properties evolve with, e.g., E or age (with 
the caveat that the evolution could be masked by effects 
arising from the geometry). Additionally, the modeling of 
pulse profiles at different wavelengths can yield constraints 
on the geometry of emission, in particular, the inclination 
angle of the magnetic axis (a) and the observer viewing 
angle (Q with respect to the spin axis. Constraints on a and 
£ can also be determined by studying the radio polarization. 
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Since the inclination of the binary system with respect to the 
line-of-sight is accurately known, £ directly measures the 
misalignment between the pulsar's spin axis and the orbital 
angular momentum, a quantity which gives important clues 
about th e formation of the b inary system (for a discussion, 
see, e.g. JFerdman et a l. 2008). 

We here report the discovery of 7-ray pulsations from the 
22.7 ms pulsar J0737-3039A, in 43 months of data recorded 
by the Fermi LAT. In Section [2] we describe the analysis of 
the LAT data and the resulting light curve and spectrum for 
PSR J0737-3039A. In Sections [XT] and £2] we model the 
offset radio and 7-ray light curves, as well as the radio polar- 
ization, under various emission models, allowing us to con- 
strain the pulsar's geometrical characteristics. We conclude 
with a brief discussion on the implication of these results. 

2. LAT ANALYSIS 

We analyzed Fermi LAT events recorded between 2008 Au- 
gust 4 and 2012 March 7, belonging to the "Source" class un- 
der the P7_V6 instrument response functions, and with zenith 
angles smaller than 100°. We rejected events recorded when 
the rocking angle of the telescope exceeded 52°, when the in- 
strument was not operating in the science observations mode 
or when the data quality flag was not set as good. Analyses of 
the LAT data were carried out using the Fermi Science ToolsQ 
(STs) v 9r28p0. The dat a were phase-folded with the Fermi 
plug-in (iRav et alj|201 lb distributed with the Tempo2 pulsar 
timing package dHobbs et al.ll2006l) . an d the ephemer i des fo r 
PSRs J0737-3039A and B published in lKramer et ail (12005) , 
Since 2008 Marc h PSR J0737-3039 B has been invisible in 
the radio domain dPerera et al . 2010) because of the preces- 
sion of its spin axis dBreton etaU2 008) causing its radio beam 
to miss the Earth, precluding timing measurements contempo- 
raneous with the Fermi mission. Additionally, changes in the 
radio p rofile of PSR J0737-3039B prevented iKramer et all 
(120061) from constructing a coherent timing model for this pul- 
sar. Our ability to accurately fold the Fermi LAT data using 
this ephemeris is therefore limited. Nevertheless, the very low 
spin-down luminosity of this pulsar makes it unlikely to be 
detected in 7 rays. 

From initial pulsation searches using standard data selec- 
tion cuts we found no evidence of 7-ray pulsations from 
PSR J0737-3039B. However, selecting events found within 
0.°5 from PSR J0737-3039A and with energies above 0.2 
GeV we o btained a value for the b in-independent //-test 
parameter (Ide Jager & Busching||2010l) of 21.3, correspond- 
ing to a pulsation significance of ~ 3.7er. The 7-ray pulse 
profile of PSR J0737-3039A comprises two peaks, which 
is reminiscent of known 7-ray pulsar light curves (e.g., 
lAbdo et all I2010al) . with no evidence for emission in the 
[0;0.2] U f0.5;0.81 phase interv al. 

As illustrated in lKeri (1201 lab and lGuillemot et al.l d2012bl) . 
weighting the 7-ray events by the probability that they origi- 
nate from the putative pulsar increases the sensitivity to faint 
pulsations. To measure these probabilities we performed a 
binned likelihood analysis of the spectra of the sources in 
the region of interest, using the PyLikelihood module of 
the STs. Our spectral mo del included the 50 2FGL catalog 
sources (iNolan et al.ll2012l) within 20° of PSR J0737-3039A. 
The contribution from the pulsar, which is not associated with 
any 2FGL source, was modeled as an exponentially cutoff 
power law (ECPL) of the form dN/dE oc £" r exp (-E/E c ), 

11 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 



TABLE 1 

Main Properties and 7-RAY Parameters of 
PSR J0737-3039A 



Parameter Value 



Rotational period, P (ms) 22.7 

Period derivative, P (lfT 18 ) 1.76 

Spin-down luminosity, E (10 33 erg s~') 5.94 

Magnetic field at the light cylinder, B LC (10 3 G) ... 4.97 

Distance, d (pc) 1150+^jj 

First peak position, <3> [ 0.43 ± 0.0 1 

First peak full width at half-maximum, FWHM[ ... 0.03 ±0.02 

Distance from the closest radio peak maximum, <5[ . 0.20 ± 0.01 

Second peak position, $2 0.93 ±0.01 

Second peak full width at half-maximum, FWHM 2 . 0- 02 ±o 02 

Distance from the closest radio peak maximum, 82 ■ 0.13 ±0.01 

7-ray peak separation, A = $2 - <J> 1 0.49 ±0.01 

Photon index, 1' < 1.3 

Cutoff energy, E c (GeV) 0.4 ± 0.4 (< 0.9) 

Photon flux, F (> 0.1 GeV) (10~ 9 photons cm~ 2 s _1 ) 6± 3 « 9) 

Energy flux, G (> 0.1 GeV) (10 -12 erg cnT 2 s" 1 ) . . 4± 1 (< 5) 

Luminosity, L 7 //n = 4irGd 2 (10 32 erg s" 1 ) 6+ 3 (< 8) 

Efficiency, 77/% = AnGd 2 /E 0. 10;*<g (< 0.13) 



NOTE. — See Section [2] for details on the measurement of these 
parameters. Quoted errors are statistical, while numbers in parenthe- 
ses represent canonical values, obtained by fixing the photon index of 
PSR J0737-3039A at 1.3 in the spectral analysis. 

where T is the photon index and E c is the cutoff energy of the 
spectrum. The source model also included contributions from 
the extragalactic diffuse emission and the residual instrumen- 
tal background, jointly modeled using the iso _p7v6source 
template, and from the Galactic diffuse emission, modeled 
with the gal_2yearp7v6_v0 map cub43 The parameters of 
PSR J0737-3039A, of the seven sources within 10° of the 
pulsar, and the normalization factors of diffuse components 
were left f ree in the fi t. A Test Statistic (TS; for a def- 
inition see No lan et all 1201 2|) value of 20.3 was found for 
PSR J0737-3039A at this stage. 

To improve the quality of the spectral results we analyzed 
the [0;0.2] U [0.5;0.8] phase interval (OFF pulse), removing 
the contribution from the pulsar. This allowed us to obtain an 
improved fit of the spectra of neighboring sources. Sources 
with TS values below 2 were excluded from the best-fit model 
obtained at this point, and the spectral parameters of sources 
beyond 3° from the double pulsar system were frozen at their 
best-fit values. We finally analyzed events in the complemen- 
tary [0.2;0.5] U [0.8; 1] phase interval (ON pulse) with the re- 
sulting source model, refitting the contributions of the pulsar 
and of the three sources within 3°, to determine the spectral 
parameters for PSR J0737-3039A listed in Table [TJ A TS 
value of 3 1 .2 is obtained for PSR J0737-3039A. Refitting the 
data with a power law shape for the pulsar, we found that the 
ECPL model is preferred by the likelihood at the 4er level. A 
cross-check of this analy sis performe d with the independent 
Pointlike analysis tool dKerr||2011bl) yielded results consis- 
tent with those listed in Table QJ Also listed in the table are 
the 7-ray luminosity L 7 = AirfaGd 2 where G is the phase- 
averaged energy flux above 0. 1 GeV and assuming a beaming 
correction factor /h = 1, and the derived efficiency 77 = L 1 /E. 
The photon flux F and energy flux G of PSR J0737-3039A 
are among the lowest values of any 7-ray pulsars detected to 
date, probably a consequence of its low E value. We also 

12 See http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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FIG. 1.— Radio and 7-ray light curves for PSR J0737-3039A. Bottom 
panel: 1.4 GHz radio profile based on observations made with the Green 
Bank Telescope. Upper panels: Fermi LAT profiles in different energy bands, 
obtained by weighting each event by its probability to have been emitted by 
PSR J0737-3039A. Photons with probabilities smaller than 0.05 were ex- 
cluded. Horizontal dash ed lines indicate the estimated background levels (see 
Guillemot et al. 1 2012b) for the determination of these lines). The 7-ray pulse 
profiles have 40 bins per rotation. Two rotations are shown for clarity. 

note that PSR J0737-3039A has the lowest value of the mag- 
netic field at the light cylinder of any recycled pulsars de- 
tected in 7 rays, .Blc ~ 4.97 x 10 3 G (see Equation 3.22 of 
lLorimer & Kramerll2005l) . 

None of the different attempts to model the spectrum of 
PSR J0737-3039A presented above yielded satisfactory mea- 
surements of its photon index, T. In all cases, a small value 
consistent with was favored. This is likely a consequence of 
the cutoff energy E c being too low to reliably determine the 
photon index V. Limits on the other spectral parameters can 
thus be placed by repeating the analysis using a value for the 
photon index which is representative of the V values observed 
for other MSPs. Table [T]lists the results obtained by fixing V 
at 1.3, this value being the average photon index measured for 
the strongest MSPs in 2PC. The spectral parameters quoted in 
parentheses represent canonical values, arising from the lack 
of constraints on V. 

The tool gtsrcprob was finally used to assign each event 
a probability that it originated from PSR J0737-3039A based 



on the fluxes and spectra obtained from the likelihood anal- 
ysis. The integrated pulse profile of PSR J0737-3039A over 
0.1 GeV and for probabilities larger than 0.05 is shown in 
Figure Q] The zero of phase is defined by the maximum of 
the first Fourier harmonic of the signal transferred back t o the 
time domain. We find a weighted H -test parameter of 41 dKerrl 
1201 latlde Jager & Busching|2010l) . corresponding to a signif- 
icance of ~5.4er. 

As can be seen from Figure Q] the 7-ray profile comprises 
two peaks, with indications for additional complexity. We fit- 
ted the profile above 0. 1 GeV with Lorentzian functions and 
found the peak positions and full widths at half-maxima for 
the first and second 7-ray peaks as listed in Table Q] Defin- 
ing 5\ and 82 as the separation between the 7-ray peaks and 
the maxima of the closest 1 .4 GHz radio peaks, respectively 
at $ ra - 0.24 and $,. 2 ~ 0.79, we find 5\ = 0.20 ±0.01 and 
62 = 0.13 ±0.01. The uncertainty on Si and 62 caused by 
the error on the dispersion delay is A(DM)/(A:/ 2 ), where 
A(DM) = 5 x 10~ 3 pc cm" 3 is the uncertainty on th e disper- 
sion measure (DM) rep orted inlKramer et al.l (120061) . k is the 
dispersion constant (see lLorimer & Kramerl2005i) and / = 1 .4 
GHz, is found to be ~ 5 x 10~ 4 in phase, and can thus be ne- 
glected. The 7-ray emission from PSR J0737-3039A is thus 
offset from the radio emission, suggesting distinct origins in 
the magnetosphere of the pulsar. Th i s also holds for the X-ray 
peaks detected by Chat teriee et al.l d2007l) in Chandra HRC 
data, which are likewise aligned with the radio peaks. 

Finally, we have searched the OFF-pulse fraction of the 
data for modulation at the orbital period, caused by the col- 
lision of the particle winds from the pulsars, and we have also 
searched the ON-pulse signal for attenuation of the emission 
from pulsar A around conjunction, caused by, e.g., photon- 
light-pseudoscalar- boson oscillation in t he magnetosphere of 
B, as proposed by iDupavs et al.l d2005j). In both cases, we 
observed only steady emission as a function of orbital phase. 

3. CONSTRAINING THE VIEWING 
GEOMETRY OF PSR J0737-3039A 

We can place constraints on the magnetic inclination a and 
viewing angle ( by modeling the radio and 7-ray light curves 
and the radio polarization. For all models we ass ume the 
vacuu m, retarded-dipole magnetic field geometry dDeutschl 
[19551) . The emission is assumed to originate in the open 
zone, determined by field lines that do not close within the 
light cylinder (at R^c = c/O = cP/2%), traced to foot points 
on the star surface at /?ns (radius of the neutron star), which 
define the polar cap. The computation and fitting procedur e 
are described in IWatters et all d2009l) . IVenter et al] d2009l) . 
iRomani & Walters! d20ldT and references therein. 

In two of the prevalent models of 7-ray emission, uni- 
formly emissive zones of the magnetosphere stand in for more 
physical mo dels. The first is the outer gap p icture (OG; 
ICheng et al.1ll986tlYadigaroglu & Romanill 19951) in which ra- 
diation is emitted between the null charge surface where 
Q ■ B = and 7?lc- The seco nd is the two-pole caustic model 
(TPC;|Dyks_&Rudak 2003), which w e take to be a geomet- 
ric re alization of the slot gap model dMuslimov & Harding! 
12001 . In the original TPC model emission extends from the 
surface to 0.75 /?lc- In both geometries the 7-ray emission is 
confined toward the edge of the open zone coinciding with the 
last open field lines, in a "gap" idealized as a region of width 
w em , interior to an accelerating layer of width w acc . In the OG 
model the 7-ray emission originates from a thin layer on the 
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inner edge of the gap, closest to the magnetic pole, whereas in 
the TPC model emission is produced throughout the gap. The 
radio emission can occupy a large fraction of the open field 
lines, albeit generally at lower altitude. 

3.1. Joint Radio and ^-Ray Light Curve Fits 

One approach is to fit the 7-ray and radio light curves 
directly, following the geometrical light curve modeling of 
I Venter et al.l ( 120091), with the fitting proced ure described in 
more detail in, e.g.. lGuillemot et alj (l2012bl) . For these com- 
putations we assume P = 25 ms and P = 1 x 10~ 18 s s -1 , and 
treat the radio pulse as a classical single-altitude, hollow- 
cone beam center ed on the magnetic axis as described by 
I Story et aTl (120071) . For this analysis we start from a grid of 
models with 1° steps in both a and f, and fit the emission zone 
accelerating and emitting gap widths, w acc and w em , with steps 
of 2.5% of the polar cap opening angle (Opc ~ (O/^ns/c) 1 ^ 2 )- 
For this computation we extend the emission to 0.95 R^c- 

The models are fit to an unweighted-counts 7-ray light 
curve for PSR J0737-3039A, constructed by selecting events 
from the data described in Section [2] found within 2° of the 
pulsar, with energies above 0.1 GeV and with probabilities 
of being associated with the pulsar larger than 0.05. In the 
case of low-statistics 7-ray light curves with sharp peaks such 
as those observed for PSR J0737-3039A, using a \~ statistic 
with binned light curves may not be optimal as this statistic 
is often insensitive to the peaks. We have thus found that it 
is preferable to use the unweighted counts, which allow the 
use of Poisson likelihood, as opposed to the weighted counts 
which require using a x 2 statistic for binned light curves. 
The radio profile is fit with a x 2 statistic, assuming the same 
relative uncertainty for each phase bin, which is combined 
with the likelihood from the 7-ray fit to determine the best-fit 
model parameters from a scan over the model phase space. 
We used 30 bins for both the 7-ray and radio light curves in 
these fits. 

One difficulty in joint fitting of the radio and 7-ray profiles 
is that the low statistical uncertainty of the radio light curve 
compared to the uncertainty of the 7-ray bins causes the fit to 
be dominated by the radio data and to ignore the constraints 
provided by the 7 rays. Because our radio model is simplistic, 
this can lead to unrealistic solutions. We have therefore in- 
vestigated different prescriptions for the "uncertainty" on the 
radio profile bins, in order to balance the radio and 7-ray con- 
tributions to the joint likelihood. The following prescription 
gives the most satisfactory results. First, we select an on-peak 
interval for the 7-ray light curve, which should dominate the 
likelihood, and define <7 gave as the average, fractional uncer- 
tainty of the 7-ray bins in this interval. We then define a radio 
uncertainty cr r = r max x er ga ve, where r max denotes the maxi- 
mum radio light curve bin. Although this prescription gives 
satisfactory results, it is arbitrary. We tested its sensitivit y 
by varying cr r by a factor of two (see lGuillemot et afll2012al) : 
the resulting best-fit geometries changed by a < 4° (12°) and 
( < 4 (9°) for the TPC (OG) model. The latter values can be 
considered systematic errors on our fit results. 

The best-fitting model light curves are shown in Figure [2] 
and the parameters are given in Table [2] The la uncertainties 
are estimated from two-dimensional likelihood profiles for a 
and ( and one-dimensional profiles for the other parameters. 
The best-fit gap widths for both models have size 0, indicat- 
ing that they are unresolved by our simulations. The simple 
geometric models we use neglect important physics, e.g., the 



effects of magnetospheric charges and currents on the gaps, 
and cannot perfectly fit the data. We use this discrepancy to 
set the scale of our parameter constraints by rescaling the like- 
lihood by its best fit value, reduced by one-half the degrees of 
freedom (dof). In the Gaussian approximation, this is equiv- 
alent to setting the best-fit reduced x 2 = 1 ■ With this normal- 
ization, the x 2 differences between the TPC and OG models 
are not significant. Also shown in Figure |2] are the closest 
approaches to the magnetic axis under the best-fit TPC and 
OG light curves. We note that our modeling finds different 
locations for the magnetic axis under the two geometries, due 
to the ambiguity in defining the first and second 7-ray peaks 
when the phase separation is close to 0.5. 

Using these geometric models, the light curves obtained 
when reflecting across a = 90° and ( = 90° are the same, 
but shifted by 0.5 in phase. We report confidence regions for 
a < 90° and ( < 90°, but note that identical, reflected regions 
exist in the other three quadrants. 

From the simu lations, we can est imate the beaming cor- 
rection factor /n (IWatters et al.ll2009l) which relates emission 
along a given line-of-sight to the total flux. The values for 
the models are given in Table [2] with estimated uncertainties 
based on the confidence regions in the a-( plane. Both mod- 
els suggest that this factor should be ~1, implying corrected 
7-ray efficiencies close to 10%. 

In addition to standard TPC and OG models, we have ex- 
plored two alternative geometries. In the first approach, a very 
low-altitude radio cone (low-altitude slot gap radio geome- 
try) was invoked in conjunction with the usual TPC model 
for the 7 rays. This was done in the context of a radio cone 
producing peaks leading the caustic 7-ray peaks. The sec- 
ond alternative, motivated by the idea that the radio may in- 
deed have a dominating leading peak with the radio profile 
lagging the 7-ray light curve, assumed conal radio and 7-ray 
geometries, with the 7-ray cone being lower than that of the 
radio. The fits from these alternative approaches were not 
satisfactory, as they were unable to reproduce the observed 
radio and 7-ray peak shapes and separations, leading us to 
abandon these scenarios. However, they point to interesting 
avenues of model refinement, e.g., investigation of cones with 
non-uniform emissivities (such as patchy or one-sided radio 
cones), or non-aligned 7-ray and radio cones. Such refine- 
ments, which are expected to lead to improved light curve fits, 
are beyond the scop e of the current pape r and will be devel- 
oped in future work (ISevffert et al.ll2013l) . 

We remark that the fitting results presented in this sec- 
tion were obtained by assuming that the observed radio emis- 
sion from PSR J0737-3039A originates from both magnetic 
poles. It has been proposed that the non-180 separation of 
the two radio peaks and the overall symmetry of the profile 
imply that the radio emission originates fr om magnetic field 
lines associated with a single pole (see, e.g., Manchester et al.l 
120051) . However. iFerdman et al.l d2013l) showed that this sce- 
nario would imply a high altitude origin for the radio emis- 
sion in the pulsar's magnetosphere. If both the radio and the 
7 rays originate from the outer magn etosphere, we wou ld ex- 
pect the pulses to be aligned in phase (IVenter et al.l2012|) . The 
two-pole origin for the radio emission therefore seems to be a 
more natural solution. 

3.2. Polarization Fitting 

When attempting to understand pulsar viewing geom- 
etry, radio polarization can provide very powerful con- 
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FIG. 2. — Top: modeled light curves and 7-ray data for PSR J0737-3039A. Bottom: 1.4 GHz radio profile and best-fitting profiles. TPC light curves are shown 
as pink lines and OG light curves are shown as green lines. The vertical dashed pink line and dash-dotted green line indicate the closest approach to the magnetic 
axis under the best-fit TPC and OG models, respectively. 



TABLE 2 

Radio and 7-RAY Light Curve Fit Parameters for PSR J0737-3039A 



Model a(°) CO w acc a (% PC ) w em (% ©pc) -ln(likelihood) dof fa 



TPC 80« 86+f 4 . . . 0.0±2.5 121.0 54 0.89^H 

OG 88+f 7 74^ O.OilO.O 0.0±2.5 123.1 53 0.95;$$ 

NOTE. — Using these geometric models the results are invariant under either a reflection across 
a = 90° orC = 90°. 

" In the TPC model the accelerating and emitting gaps are the same. 



straints. Although recycled pulsars are notoriously difficult to 
model, PSR J0737-3039A shows substantial linear polariza- 
tion structure, and thus offers good prospects of constraining 
the orientation. In our polarization modeling we follow the 
conventional assumption that the electric vector position an- 
gle (P.A.) follows the projection onto the plane of the sky of 
the magnetic field line at the emission point. The P.A. sweep 
as one moves past the radio pole(s) thus probes the viewing 
geometry, with particularly strong constraints on j3 = (-a. 
Our modeling extends beyond the point dipole rotating vector 
model (RVM; iRadhakrishnan & Cookdfl969l) . to follow the 
distortion of the magnetic field as t he emission zone move s 
to a non-negligible fra ction of R] p. (Craig & Romanill20l2h . 
In addition, we follow iKarastergiou (120091) in treating "or- 
thogonal mode jumps" for which the P.A. shifts by ±90°, 
and optionally account for the effects of interstell ar scattering 
neglig ible for PSR J0737-3039A). For details see lCraig et all 
200) . 

Figure [3] shows a polarimetric profile for 
PSR J0737-3039A measured at 1.4 GHz with the Parkes 
radio telescope (R. N. Manchester, private communication). 
The Stokes parameters and Gaussian decomposition of the 
linear intensity components are displayed in the top panel, 



TABLE 3 

Polarization Fit Parameters for PSR J0737-3039A 



a f (°) 


C/(°) 


n/Ru: 




dof 


At/> 


98.8+f 5 


95.8+' 3 3 2 


u - -0.01 


0.11*$ 48 


35 


n 44^+0.008 



NOTE. — Errors are the extrema of the lcr contours in the full mul- 
tidimensional parameter space. T| is the emission altitude of the central 
component of PI and rj the altitude of the central component of P2. Aifi 
marks the offset of the closest approach of the surface magnetic axis (in 
PI) from the total intensity peak in P2; this is phase 0.243 in Figure[T] 



while the lower panel shows the P.A. data. In this section we 
refer to the brightest radio peak as P2. Both it and PI are 
multi-component. In polarization fitting, these components 
are interpreted as different emission zones with different 
linear and circular polarization fractions and, possibly, 
different orthogonal mode states and emission altitudes. 
Progressively across the PI pulse we see a component 
with strong polarization at nearly constant P.A., a largely 
unpolarized peak, a component with rapid linear polarization 
sweep after the maximum and a separate, weakly linear 
polarized component on the flat pulse tail. For P2 the pattern 
appears reversed, with a weak linear tail leading the pulse, 
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a linear component with rapid P.A. sweep, an unpolarized 
peak and then a strongly linearly polarized component with 
little P.A. sweep. The flat P.A. components at the front of 
PI and the back of P2 are not naturally produced by any 
component locked to a local magnetic field. We suspect that 
these components are controlled by magnetospheric plasma, 
either in fixing the P.A. directly or in strongly distorting the 
field at very high altitude. Here we fit the rapid P.A. sweep in 
the central linearly polarized component of each pulse. The 
parameters are a, (, </>o (the phase of the total intensity peak 
with respect to the surface dipole magnetic axis), 'I'o, the P.A. 
offset and the two emission heights r\ and r2 in units of R^c- 
The fit values, the \ 2 an d the full projected (multi-parameter) 
uncertainty ranges on the fit parameters are shown in Table [3] 
Figure |4] shows a projection of the \ 2 surface of the fit onto 
the a-( plane. 

The best-fit P.A. models for the two poles are displayed in 
the lower panel of Figure [3] Note that the weak trailing linear 
component in PI is separated from the central component by 
a zero in the linear intensity L, and so is very likely an or- 
thogonal mode jump. The jumped P.A. points cluster around 
the model prediction. We do not attempt to model the total 
intensity light curve or the P.A. behavior of the flat P.A. com- 
ponents leading PI and trailing P2, since these do not follow 
simple magnetospheric models. However, we can make some 
comments about the origin of this emission. First, the circles 
mark the phase ranges of the open zone at the altitude fit for 
the two central components. These components occupy the 
second half of the open zone at this altitude, terminating at the 
closed zone boundary, an excellent consistency check. How- 
ever, the wings of the pulses, especially at large <fi, lie beyond 
the circles and hence must either come from higher altitude 
or else arise from the classical closed zone. Our picture of 
each pole is thus a cut through a hollow cone with the cen- 
tral component arising from low altitude and higher altitude 
components from the pulse wings. Since PI has the smallest 
\/3\ = 3°, it is viewed at lower altitude than P2 with /3 = 14.°6. 

These fits also determine the phase of closest approach of 
the surface dipole axis. This can be compared with predic- 
tions of the ou ter magnetosphere 7-ray m odels at the fit a and 
£. We follow iRomani & Wattersl ( 120101) in examining the a 
and £ for the best fits to the 7-ray light curve and then check- 
ing against polarization constraints. For the OG model the 
best light curve matches are found very close to the angles fit 
to the polarization data. For the TPC model (here limited to 
0.75/?lc), the best fits also occur near (90°, 90°). These are 
somewhat worse than the best OG values, since they include 
too much unpulsed emission at all phases, but they are quite 
acceptable. In this regard the polarization fitting agrees well 
with the fits of Section [3~Tl when transformed to the a > 90° 
and ( > 90° region. However, the very accurate phase mea- 
surements from the P.A. sweep highlight a tension between 
the two bands: the radio-determined pole is always ~ 0.1 in 
phase before the best-fit axis for the 7 pulse. This tension is 
also evident in the fits to the cruder pulse profiles of Section 
13.11 However, non-vacuu m effects provide a plausible expla- 
nation for this offset. In Kalapotharak os et al.l (120121) it was 
noted that the primary effect on the high altitude pulse of in- 
creasing the magnetosphere conductivity was a growing lag 
in phase with respect to the surface dipole. Lags of 54> ~ 0.1 
were seen for high conductivity magnetospheres. If we dis- 
connect the 7-ray and radio zero-phase definitions in the joint 
light curve fits, letting 6 be larger by as much as 0. 1 in phase, 



the best-fit a and £ values changed by < 4°, generally making 
the agreement between the two methods better. 

In Figure [5] we show the predicted OG and TPC 7-ray light 
curves at the pulsar geometry determined from the polariza- 
tion study, superimposed on the radio and the 7-ray data, and 
using the phase of closest approach to the surface dipole axis 
as defined by the polarization fits. Both modeled light curves 
are acceptable, with the OG one being somewhat better (re- 
duced x 2 °f 1-45 compared to 1.56 for the TPC model). For 
these a and £ angles we find beaming correction factors /n of 
0.92 for the OG model and 1.15 for the TPC model, similar 
to the ones given in Tableland therefore also implying 7-ray 
efficiencies of ~ 10%. 

4. DISCUSSION 

4. 1 . The Low Spectral Cutoff Energy and Magnetic Field at 
the Light Cylinder ofPSR J0737-3039A 

In all outer-magnetospheric models, we expect the acc eler- 
ating electric field E\\ to scale as (e.g.. lAbdo et al.ll2010bl) : 

E,| (xC(r)B LC wL, (1) 

with C(r) some function of radius, and w acc the accelerating 
gap width. On the other hand, if leptons are accelerated in 
the radiation-reaction-limited regime, i.e., when the gain in 
energy due to acceleration by the electric field is balanced by 
curvature radiation losses, we expect that the spectral cutoff 
energy should scale as: 

E c KEfp\I^C{r?I^McP l J\ (2) 

with p c the radius of curvature (at the position where the 
photons of maximal energy are produced). If we take C(r), p c 
and w acc constant for all MSPs (e.g., independent of P and P), 
this leads us to expect: 

E c ocB^. (3) 

One does see a trend when plotting published values of 
E c versus Z?lc for all known gamma-ray MSPs, but it is 
much weaker than expected, and the index is much smaller 
than 0.75, implying that our assumption of constant C(r), 
p c and Wacc across the MSP population is too simplistic. 
PSR J0737-3039A is a mildly recycled MSP, having a rel- 
atively long P and hence the lowest Z?lc oc p~ 5 l 1 p x l 1 of all 
currently known gamma-ray MSPs. Assuming that the above 
empirical trend may be extrapolated down to low Z?lc values, 
one would expect a relatively low cutoff energy. This is in- 
deed seen, although PSR J0737-3039A's cutoff energy falls 
significantly below even the extrapolated trend. 

Apart from invoking the weak empirical trend discussed 
above, one may more generally observe that the relatively 
long period of PSR J0737-3039A implies a lower accelera- 
tion potential and hence a lower spectral cutoff than for typi- 
cal MSPs, as is observed from the spectral analysis. 

4.2. Constraints on the System's Formation History 

From the modeling of the radio and 7-ray light curves of 
PSR J0737-3039A under the OG and the TPC geometries, 
and the modeling of radio polarization data (see Section [3]), 
we constrained the magnetic inclination angle a and the view- 
ing angle £ to be close to 90°. The best-fit a and £ angles 
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FIG. 3 . — Polarimetric profile for PSR J0737— 3039A as observed at 1 .4 GHz with the Parkes radio telescope (R. N. Manchester, private communication). Top: 
Stokes parameter curves (black: I, red: L, blue: V) and Gaussian decomposition of the linear intensity components (dotted curve is all linear, the thick green 
curve gives the rapid sweep central component fit here). Bottom: position angle data (blue: central components, yellow: other PA. values, green: PI tail with an 
orthogonal mode jump). The smooth curves give the best fit model for the two poles while the red circles denote the boundaries of the open zone at the emission 
altitude. Arrows denote the phase of the closest approach of the magnetic axes to the Earth line-of- sight. The smallest impact parameter /3 is for PI at phase 
<Pb = (phase of the magnetic dipole). 




FIG. 4. — Projection of the \ 2 surface of the RVM fit to the radio polar- 
ization data (central components) for PSR J0737— 3039A onto the plane. 
The best fit orientation angles (a, Q are indicated by an arrow, while the 
contours are shown at the 1<j, 2cr and 3<r levels. 



obtained from the polarization study (see Table [3} agree well 
with those obtained from the light curve modeling (Table 0. 
The conclusio n that a and C are clo se to 90° is in line with 
the results of iFerdman et al.l (12013b , based on the analysis 
of six years of radio observations of PSR J0737-3039A, re- 
vealing no significant variations of its radio profile with time. 
IFerdman et all (120131) find a and < values of 90?2+\ 6 6 3 2 and 
90.°8 

-a46' respectively. These values were calculated assum- 



ing that radio emission is seen from both magnetic poles, that 
the orbital inclination is 88. °7, and averaging the results from 
all radio pulse heights used (30% to 50%) and all observation 
epochs. 

By virtue of the relationship between £ and the misalign- 
ment angle between the spin axis and t he orbital angular 
mom entum, 8$o (cf. Equation (3.36a) of iDamour & Tavloil 
fl99l . our values for £ close to 90° imply small 6 so values, 
unless the pulsa r is at a special preces sion phase. Neverthe- 
less, the fact that lFerdman et al] (120 131) find a 95% confidence 
upper limit on 6$o of ~ 3.°2 from their observations taken 
between 2005 June and 2011 June indicates that the pulsar 
was very unlikely to be at a special precession phase dur- 
ing the Fermi observation. Our constraints on ( therefore 
provide an independent confirmation that 5 so is close to 0°. 
Since the spin axis of pulsar A and the angular momentum 
of the binary are likely to have become aligned through the 
accretion of matter from pulsar B's progenitor star by pulsar 
A, a value close to 0° for the misalignment angle Sso sup- 
ports the scenario under which a small kick was imparted 
to the system by the supernova explosion of pulsar B's pro- 
genitor star. An electron-capture supern ova, resulting from a 
collap sed O-Ne-Mg core as proposed bv lPodsiadlowski et al.l 
d2005|). could provide su ch a formation scenario for pulsar B 
(see lFerdman et al.l2013l and references therein for a detailed 
discussion). 

Because PSR J0737-3039A is a faint 7-ray emitter, we 
were not able to investigate any evolution of the 7-ray profile 
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Fig. 5. — Light curves for PSR J0737— 3039A at the a and (, angles determined from the radio polarization study, for the OG model (solid blue line) and the 
TPC model (dashed green line). The observed radio and 7-ray profiles are shown as a solid red line and as a shaded gray histogram, respectively. See Figure[3] 
for the definition of the magnetic axis phase <f>B- The blue arrow denotes the location of the magnetic axis under the OG and TPC geometries. 



as a function of time caused by a variation of the line-of-sight 
angle, £, in the four years of LAT data considered here. Since 
the misalignment angle 6$o has been shown to be close to 0°, 
changes in £ are likely small, but are still expected. An ex- 
citing prospect, when significantly more Fermi LAT data are 
available, will be to monitor the evolution of the 7-ray profile 
of PSR J0737-3039A. To date, no 7-ray pulsar light curves 
have even been observed to vary with time. 
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